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Abstract 

Many experiments on high-temperature superconductors have shown paramagnetic behavior 
when the sample is field cooled. The paramagnetism was first attributed to a ti-wave order pa- 
rameter creating vr-j unctions in the samples. However, the same effect was later discovered in 
traditional low-temperature superconductors and conventional Josephson-j unction arrays which 
are s-wave. By simulating both conventional and mixed vr/conventional Josephson-j unction arrays 
we determine that differences exist which may be sufficient to clearly identify the presence (or 
absence) of vr-j unctions. In particular the vr-j unctions cause a symmetry breaking providing a 
measurable signature of their presence in sample. 

PACS numbers: 74.50+r, 74.72-h, 75.20 -g 
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In recent years a significant discussion in the superconductivity field was devoted to the 
paramagnetic response in field-cooled samples, the so-called paramagnetic Meissner effect 
(PME), |T]|. PME was first observed in high-T^ ceramic materials [|T], 0, ^ and then also 
in low Tc samples A mechanism to explain the effect was introduced by Sigrist 

and Rice where PME was due to the presence of vr-junctions in the multiple-connected 
network of junctions formed by the superconducting grains in the high-Tc ceramic, vr- 
junctions are Josephson junctions formed between superconductors with unconventional 
pairing which cause a vr shift in the phase-current relation |^. Unconventional pairing, 
such as a d-wave order parameter, was introduced to explain the properties of high-Tc 
ceramic materials 0. PME was taken as strong evidence for a d-wave order parameter 



T0| . Multiple-connectiveness and vr-junctions were soon shown to give a paramagnetic 



response in simulated Josephson junction arrays by Dominguez, Jagla and Balseiro [11 



The discovery of PME in conventional low-Tc (LTC) samples shows that PME cannot 
always be attributed to vr-junction because they are present only in d-wave superconductors 
(or in unconventional superconductor- ferromagnet systems H). New explanations such as a 



giant flux state [12|, flux compression [13 and surface states \T^] have been introduced. The 
common feature was that PME in LTC samples is described as non-equilibrium phenomenon. 

Very recently a new experiment was devised to test the connection between multi- 
connectiveness and PME in conventional systems. A square array of LTC junctions was 
fleld cooled and shown to be predominantly paramagnetic [T3|; this paper also proposed a 
qualitative explanation for the effect. A subsequent analysis performed by numerical simu- 
lations of the arrays conflrmed the presence of PME in LTC arrays [jTB]. On this basis the 
essential ingredients of PME appears to be multiple-connectiveness rather than the presence 
of vr-junction or non-equilibrium effects. 

It is clear that simulated Josephson-j unction arrays both with[jlT| and without [[1^ vr- 
junctions show PME. Here we investigate by means of numerical simulations the differences 
between the two cases and the true role of the vr-junctions. We also suggest an experimental 
signature of the presence of vr-junctions in fleld-cooled array samples. 

The array is described by means of a full mutual inductance model similar to that used 



in ref. [|T6|. The main difference is the presence of some vr-junction, i.e., Josephson junctions 
where supercurrent is proportioinal to sin(y9 + vr) instead of simp (cf. Fig. 1). The dynamics 



of an NxN array are described by the following equations [16, Qq, Q 
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^ (1 + 5) sm{0 + ^) + J^l^ + 1$ = kL~^m. (1) 

iTT V PC 

Here time is normalized to a cell frequency (ci^"^ = L'C), with L' the self inductance of 
the elementary loops formed of the array C the junction capacitance. jSc is the Stewart- 
McCumber parameter, f3c = 27r/oC/$oG'^ where Jq is the (mean) Josephson critical current 
for the junctions in the array, $o is the flux quantum and G is the junction conductance. (3l 
is equal to 2ttL'Iq/^q. m represents the normalized loop magnetization (in unit of $o) which 
depends on frustation /, phase vector (p and the initial distribution of flux quanta in the 
array loops n, here choosen to be a random integer vector |T^. K and are the matrixes 
carrying the mutual inductance coupling. The vector ip has values of or vr, providing a 
TT-shift to a subset of the junctions in the array. The concentration c of vr-j unction |Tl|] is the 
ratio between the number of vr-junctions and the total number of junctions. The vector 5 is 
a Gaussian variable with zero mean which we use to introduce disorder in the distribution 
of critical currents. 

Here we suppose that in the array the contribution of the screening currents is not neg- 
ligible, i.e., /3l > 1 (cf. [^). In particular we choose (3l = 30 and (3c = 63 (cf. refs. 
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In both experiments [15] and previous simulations [16| it was shown that multiple- 
connect iveness effects are dominant for large field, i.e., for large values of frustation /. 
To understand the effect of the vr-shift we studied smaller values of frustation < / < 1. In 
this region the behavior of an array without disorder and with all normal or all vr-junctions 



roughly follows the single loop behavior |15, 16, T7]which is described by 

J = — sm — n (2) 

$0 27r ^2 4 ^ ' 

where $ is the total flux throught the loop. Solution of Eq. (2) gives the possible states 
for the total flux $ (or the current) in the single loop. For if) = {] Eq. (2) describe a single 
loop made by four conventional Josephson junction. There are four independent solutions 
which are obtained varying the integer n from zero to three. The (stable) lowest Gibbs 
energy state is diamagnetic for < / < 1/2, i.e., the loop current is negative, generating 
a magnetic moment opposite to the external fleld, and paramagnetic for 1/2 < / < 1, i.e., 



the loop current is positive |15, 16]. On the other hand li ip = the vr-shift will reverse 
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the behavior: for the vr-loop (a loop with an odd number of vr-junctions) the lowest Gibbs 
energy state is paramagnetic for < / < 1/2 and diamagnetic for 1/2 < / < 1. 



In Ref. |T6| we showed that conventional-array behavior can be qualitatively understood 
in terms of a model based on the single-loop behavior. We found that each loop in the array 
had a magnetization that was close to one of the lowest energy solutions of Eq. (2) for the 
isolated loop, with the number of loops of each type determining the overall magnetization. 
This simple picture also works for arrays with both conventional and vr-junctions, with the 
number of possible solutions being increased because that Eq. (2) has additional solutions 
for the vr-loops. 

In order to compare conventional arrays with mixed conventional/vr arrays, we first sim- 
ulate the conventional arrays using the method of Ref. |jl6[ . In Fig. 2a the magnetization of 
a 20 X 20 conventional array with / = 0.2 is shown. There are two single-loop magnetization 
states with a predominancy of diamagnetic states (light gray) over which a set of paramag- 
netic loops (dark gray) are randomly distributed. The histogram of loop magnetizations is 
shown in Fig. 3a and shows clearly only two peaks corresponding to two loop magnetiza- 
tions. The mean magnetization is negative and small as expected for such small values of 
frustation accordingly to the Eq. (2) for the single loop; the exact value is m = —0.0027. 
As said above the values of the single peak averages can be obtained from Eq. (2); the small 
spread of the peaks is due to mutual inductance coupling. 

In Fig. 2b a similar 20 x 20 array with / = 0.2 is shown in which 380 vr-junctions have 
been introduced at random. This is equivalent to a concentration of vr-junctions of c = 0.45. 
The main difference with the case of coventional array is the presence of four magnetization 
states, two for each type of loop. We note that the new vr-loop states are predominantly 
paramagnetic according to the symmetry shift discussed above. This, together with the 
lower number of conventional loops, is causing the reversal of the mean magnetization (here 
being m = +0.0018). This is the way in which a mixed array becomes paramagnetic for 
low values of frustation. Lower concentrations will make the array diamagnetic (i.e. the 
dependence on concentration is similar to that reported in Fig. 3b of Ref. |jTl[ for a 3D 
system). For other values of / we note that the single loop symmetry is broken, i.e., for 
/ < 1/2 the array can be diamagnetic or paramagnetic depending on the distribution and 
number of vr-loops. 

From the preceding discussion, it is clear that there are at best only subtle differences 
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between a conventional array and an array containing a mixture of conventional and vr- 
junctions. In the remainder of the paper we will show that when / = 1/2 these differences are 
the most pronounced, and therefore / = 1/2 is probably the best place to look experimentally 
in searching for evidence of vr-junctions in non zero magnetic field. When / = 1/2 the 
problem is simplified by the fact that the four lowest states shrink to three because the 
vr-loop magnetization becomes zero. So one is confronted with only three peaks which are 
symmetrically distributed around the zero magnetization. In Fig. 4a this case is shown in 
a magnetization histogram. The central peak is the zero magnetization of vr-loops and two 
symmetric lateral peaks are due to coventional loops. Symmetrical histogram structure is 
an optimal "signature", it is easily recognized with SSM and statistical filtering or image 



filtering techniques can be applied knowing the true distribution |21 



The central peak intensity is related to presence of vr-loops. At c = 0.3 (250 7r-junctions) 
it becomes thinner preserving the same height. Then it decreases by about 10% at c = 0.15 
(126 TT-j unctions). Finally at c = 0.05 (42 vr-junctions) it decreases of about 60%. 

In Fig. 4a a weak Gaussian disorder with a standard deviation of cr = 20% was introduced 



in the critical currents of the arrays similarly to Ref.[20|. This was added to shown that 
weak disorder does not change the three-peak structure of the magnetization histogram. 
For high values of the parameter (Sl the nature of the magnetization states is discrete and 
this is preserved for weak disorder in the currents. Peak structure exists also with stronger 
disorder in the critical currents. In Fig. 4b we set a = 80% |^ we see peaks enlarging until 
the tails of the single three distributions overlaps. It is observed that the remaining discrete 
structure is due again to the large jSi loops. Also we have simulated the array with the same 
parameters of Fig. 4b using a mean value of 1/3 for the vr-j unction critical current. The 
result is again similar to Fig. 4b but that central peak is lowered by roughly 13% and the 
histogram appears more enlarged. 

In conclusion we would stress again how both conventional and mixed conventional/vr 
arrays show a similar response when field-cooled. This response is related to the single loop 
solutions of Eq. (2). In low field / < 0.5 paramagnetism in mixed conventional/vr arrays is 
the result of the paramagnetic behavior of 7r-loops. To quantify this effect for a disorderd 
system is difficult because it depends on distribution and number of vr-loops. On the other 
hand a noteworthy point is / = 1/2. Here the vr-loops sets to zero current so it is possible 
from a measurement of the magnetization histogram to trace back the presence of vr-loops in 
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the sample. For High-Tc materials this experiment can be an important complement of the 



search for spontaneous currents in zero field Correlation between / = and / = 1/2 
would be of great interest permitting to estimate the vr-loops content of the sample. 
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Figure Captions 



Fig. 1 Two-dimensional Josephson junction array with a random distribution of tt- 
junctions. The gray loops represent the 7r-loops, i.e., loops with an odd number of tt- 
junctions. The figure is just the top-left corner of the mixed array in Fig. 2b. 

Fig. 2 Simulated magnetization of a 20 x 20 array at / = 0.2 with = 30 and 13c = 63: 
a) without TT-junctions, light gray represents diamagnetic loops, gray paramagnetic ones; b) 
with TT-junctions at a concentration of c = 0.45 white and light gray represents diamagnetic 
loops, gray and dark gray paramagnetic ones. 

Fig. 3 Magnetization histogram of a 20 x 20 plane array at / = 0.2 with I3l — 30 
and (5c — 63: a) distribution without TT-junctions; b) distribution with TT-junctions at a 
concentration of c = 0.45. 

Fig. 4 Magnetization histogram of a 20 x 20 plane array at / = 0.5 with jS^ = 30, 
Pc = 63 and TT-junctions distribution at c = 0.45: a) effect of a weak disorder in the current 
distribution {a — 20%); b) the same for a stronger disorder (cr = 80%). 
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Fig.2 De Leo & Rotoli 
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Fig.3 De Leo & Rotoli 
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Fig.4 De Leo & Rotoli 



